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The theoretical concepts for the renormalization of self-consistent Dyson resummations, 
dcviced in the first paper of this series, are applied to first example cases for the (/)^-theory. 
' Besides the tadpole (Hartree) approximation as a novel part the numerical solutions are 

, presented which includes the sunset self-energy diagram into the self-consistent scheme based 

' on the ^-derivable approximation or 2PI effective action concept. 
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INTRODUCTION 



I In the first paper of this series Q (in the following referred as I) we have derived the theoretical 

' concepts for the renormalization of Dyson equation based resummation schemes at finite temper- 

■ ature. It could be shown that such truncated self-consistent approximations can be renormalized 

with counter term structures solely defined on the vacuum level of the considered approach, if two 
conditions are met: a) the underlying exact theory has to be renormalizable and b) the appro- 
ximation scheme has to be based on Baym's ^>-derivable concept |3|, i.e., on a two-particle 
, irreducible (2PI) effective action principle. Thus the self-consistent self-energies S are generated 

\ from a truncated set of 2PI closed diagrams (called ^) of the underlying full theory through 



where is a functional of the self-consistent propagator G{p) in momentum representation (as 
^ . in I an asterix behind an equation number implies that the corresponding relations are only valid 

^ I at the regularized level, while all other equations are valid also for the renormalized quantities). 

To repeat, the main issue is not to render any divergent loops finite, this has been pursued many 
times. The aim is to deploy the counter term structure such that it is entirely determined at the 
vacuum level of the self-consistent scheme. Besides the explicitely visible divergences this implies 
to resolve also the nested and overlapping vacuum divergences hidden in the self-consistent matter 
parts of the propagator. 

In I we have given the proof that for the self-energy this can be done for any <I>-derivable 
approximation, provided the underlying quantum field theory is renormalizable in the usual sense. 
Thereby the counter-term structure results from closed equations on the vacuum level, implicitly 
generating a particular though infinite subset of counter terms. The complexity of the ensuing 
equation is similar to that for the self-energy. 

Thereby it is of particular importance that the entire counter term structure is consistently 
constructed solely and only within the effective action defined by the chosen approximation to <I>. 
This implies that, e.g., the so obtained counter-term scheme and thus the behavior of the running 
coupling constant clearly deviates from the true one at orders beyond those included in If 
these rules are not watched one may face uncompensated divergences as recently encountered in 
the detailed work of Braaten and Petitgirard They tackled the same problem, however using 
different techniques in terms of a restricted ansatz for the propagator. Uncompensated singularities 
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arose from the fact that the used /5- function was not the one pertaining to the self-consistent scheme 
but taken from the fuh /3- function calculated up to the fifth loop order. Furthermore it is yet not 
transparent to us, how in their approach the hidden and nested vacuum sub divergences can be 
resolved such that one arrives at a counter term structure solely defined on the self-consistent 
vacuum level (T = 0). In our case the latter leads to a to be renormalized Bethe-Salpeter equation 
for the vacuum four-point function consistent with the chosen approximation level for the $- 
functional. 

In I we also have shown that the generating functional r[(/3, G] and thus the thermodynamic 
potential can be rendered finite with counter terms solely defined on the vacuum level of the 
self-consistent scheme. Again it is crucial to expand the functional around the solutions of the 
corresponding self-consistent approximation for the vacuum. It is not possible to render arbitrary 
parts of r finite for themselves since it is important to use the stationarity of F for the solution of 
the vacuum equations of motion. In particular terms linear in the matter parts of the propagator 
are singular and only drop out if the vacuum part is solved self-consistently! 

For self-consistent Dyson resummation schemes it is well known that they may violate 
symmetry properties such as crossing symmetry, masslessness of Nambu-Goldstone modes, conser- 
vation laws, etc. at a level beyond that of one-point functions, i.e., on the correlator level at orders 
beyond those included in $. We devote a forthcoming paper 0] in this series to the issue of how 
to extend the scheme such that symmetry violations are indeed cured. 

In this paper we give first numerical solutions for the leading and the next-to-leading order of 
the self-consistent Dyson equations for the cp^-theory, which is defined by the Lagrangian 

^ = l{d,cP){d^^)-^^'-^cp\ (2) 
In the symmetric phase, {(j)) = 0, the first diagrams contributing to $ and S are 

^= OO +1 +•••' (3*) 




(4*) 

where all lines represent self-consistent propagators. The leading order diagram gives rise to the 
tadpole (Hartree) approximation for the self-energy, which, frequently considered in the literature 
(see for instance [^, |8|), leads to the standard gap equation. Here it is given as the most simple 
example for our general renormalization scheme. The next-to-leading order includes the sunset 
diagram for the self-energy. While in the tadpole case the self-energy is real and constant, with 
the sunset term the self-energy becomes momentum dependent and complex. Thus the particles 
acquire a finite spectral width due to collisions with the surrounding matter. With the sunset 
diagram one enters a new stage of sophistication, both, as far as the counter term structure is 
concerned, and with respect to the numerical solution of the resulting self-consistent equations. 

First numerical investigations of the self-energy are presented which result from the two explic- 
itly given diagrams in (13). We could improve the numerical accuracy immensely compared to the 
status given in where we restricted ourselves to the computation of the imaginary parts of the 
self-energy to avoid problems with renormalization. The paper is organized as follows: In sect. || 
we briefly summarize the results of I as far as we need the formulae for the numerical calculations. 



In sect. |ll| we solve the tadpole approximation (gap-equation) and in sect. IV the next-to- leading 
order approximation including the sunset diagram. 

Throughout the paper we work in momentum-space representation within the real-time field 
theory formalism for the self-energies and propagators. Thus, if not stated otherwise, self-energies 
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and Green's functions are contour matrix valued in the sense of the Schwinger-Keldysh real-time 



contour |10| , 11, [T^. There vertices belonging to time arguments on the time and anti-time ordered 
branches are labeled by — and -|- superscripts, respectively. The integration sense is accounted 
for by assigning the values itiA to the bare vertex depending on its placement on the {— } or {-|-} 
branches of the contour. 



II. THE RENORMALIZED EQUATIONS OF MOTION 

In I we have shown that arbitrary self-consistent approximations given by the ^-derivable con- 
cept (0) can be renormalized with counter terms defined on the vacuum level. Thereby the 
renormalized self-energy as well as the thermodynamical potential are finite and consistent to one 
another. 

For the renormalization procedure the self-energy has to be split into three parts, namely the 
pure vacuum part (which is of divergence degree 6{T,^^^'^^) = 2), the part T,^^^ with divergence 
degree (5(S('^)) = 0, containing explicit and hidden vacuum sub-divergences, and the rest S^''^ with 
6{T,^^^) = —2, which contains only explicit vacuum sub-divergences with corresponding counter- 
terms 

S = + + (5) 

The self-consistent Green's function G follows from the Dyson equation in contour matrix form 

A-^G=1 + J:G. (6) 

Here A is the free Green's function on the contour. In view of (|5|) the self-energy within the Green's 
function needs two subtractions to render its loops finite: 



iG = iG^™") + iG'(™^)s(°''^'^)G'(™^) + iG^'') = + ^^^^^^1 + ■ 



(7) 



j(^(mattcr) 



The subtraction only affects contour diagonal pieces of the propagator, since loops containing mixed 
contour vertices are finite for themselves. Thus 

1^ iP)-[ iG++(-^)(p) J ' iS++(0)(p)J' 

and likewise S'^™'^) are diagonal. While G^™'^^ and S'^™'^) are given by the (anti-)time ordered self- 
consistent expressions given by the ^-derivable scheme on the vacuum level, T,^^^ accounts for the 
self-energy parts linear in Qi'^^^^^^) _ The three r.h.s. terms of (|7|) are of divergence degree —2, —4, 
and —6, respectively. It is important to notice that only the full self-consistent propagator and self- 
energy, G and S, obey the equilibrium conditions (KMS) among their four contour components, cf. 
Eqs. (A. 14) to (A. 18) in I. The components of subtracted pieces like G^™^"''''), S^^^ or S^'') obey no 
direct equilibrium relations as they result from differences of finite temperature with vacuum terms. 
Here we see an advantage in the real-time formalism which naturally permits this decomposition 
into vacuum and matter pieces of all dynamical quantities. 

The contour diagonal parts of T,^^^ and T,^^^ contribute only at finite temperature. Due to 
Lorentz invariance the vacuum self-energy and the Green's function are functions of s = p^. The 
vacuum self-energy is renormalized according to the rules of the on-shell renormalization scheme, 
i.e., 

S(™=)(s = m2) = 0, 9,S(™'=)(s = ?n2) = 0, (9) 
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which defines m to be the physical mass of the particles and normalizes the wave functions such 
that the residuum of the Green's function at s = is equal to 1. As we shall see in the following 
the renormalization procedure used to calculate the vacuum parts, and subsequently needed for 
the temperature dependent pieces of the self-energy, allows to chose any renormalization scheme. 
Especially it is possible to use "mass independent" minimal subtraction schemes (MS [^] or MS 
Ijl^ of dimensional regularization) so that the renormalization program shown in I is applicable 
also for the massless case. For the renormalization of the self-energy and the 2PI- functional T[G] it 
is crucial to split all quantities into vacuum (T = 0) and finite temperature parts by the procedure 
given above. 

As shown in I this procedure allows to define the generating functional in its dependence on 
the in-matter parts of the propagator with the self-consistent solution given by the stationary 
point. Simultaneously this leads to the renormalized thermodynamical potential. All divergences 
indeed compensate if and only if the vacuum level is solved and renormalized entirely within the 
same approximation level as used for the finite temperature parts. This fact together with the 
stationarity conditions of the generating functional at T = for the vacuum equations of motion 
avoids the singularity problems encountered in ||5|. 

In S^''-* the in-matter part of the propagator is involved in logarithmically divergent loops which 
imply hidden divergences. As shown in I these can be resolved with the following result 

^^'Hp) = I ^{[rW(/,p)-r(4'™^)(/,o)]G(--"-)(/) + A(o,/)G«(/)}. (10) 

Here in consistency with the Dyson equation the four-point function T^'^\p,q) is generated also 
from $ through 



The counter term T^^'™'^) contains only the contour diagonal parts (all vertices placed on one side of 
the real-time contour) of T^^^ . This function is 2PI in the channel p ^ q, i.e., one cannot disconnect 
the diagrams by cutting two lines and separating the lines carrying the momenta p and q. It 
contains only vacuum pieces and can straight forwardly be renor malized with r(^'™'=)(0, 0) = ±A/2 
(denoted by the superscript (ren)). At the same time this four-point function defines the kernel 
for the vacuum s-channel Bethe-Salpeter equation defining the 2PR four-point function A{p,q). 
In (Il0|) we only need the half sided A(0,p), which, as one of the crucial points of the Bogoliubov- 



Parasiuk-Hepp-Zimmermann (BPHZ) renormalization procedure |15, [T^, is given by the finite 
equation 

A(o,p) = ±^ + r(4'-<=)(o,p) + i I ^A(o,0[G(0]'[r(4'^-'=)(/,p)-r(4'™=)(/,o)], (12) 

again involving only contour diagonal terms. Here we inserted the renormalization condition 
A(0,0) = ±|. 

The expression ( [lO|) for S'^'^^ void of the A(0, /) part would correspond to a naive subtraction, 
which indeed is finite, however it implies T dependent counter terms. As explained in I, only 
the entire expression (|^) guarantees a counterterm structure solely defined on the self-consistent 
vacuum level. 

The perturbative view 

The above procedure defines tiie self-consistent renormalization scheme in terms of the corresponding self- 
consistent vacuum propagator the vacuum four-point function A(0, 1), and the self-consistent in-matter pieces 
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of the propagator (;;(™atter) As an illustration we briefly quote the corresponding perturbative (fPI) view of 

the scheme, expanding the full self-energy in terms of the vacuum propagator in contour matrix notation satisfying 
the equilibrium relations (KMS) at the given temperature T 



.g{vac;T). . _ / iG '^^^P) ] + o'^^''^'> ( \vn\ v) f "^^'^^ 0{~Po) + n{po) 



(13) 



+ 



with the vacuum spectral function and thermal Bose-Einstein weight defined by 

^ -2lmG^r\p) = -2sign(po)ImG— 

1 (14) 
- exp(/5|po|)-l- 

As discussed in I, loops involving the O-functions or the thermal weights summarized by the dotted line are always 
finite. 

For the self-consistent scheme of tadpole and sunset the following types of IPI diagrams are resummed to the 
full self-consistent self-energy: 




+ ) +••• + —♦ +■■■ + — 4 ^ -f (15*) 



y^(vac) 



MO) 



-1 (matter) 



Here the diagrams for generated by iterative insertions of the diagrams of Y,^^'^'^') ^ replacing one or more 

plain lines by dotted ones with the constraint that due to the self-consistence at the vacuum level all pure vacuum 
self-energy insertions are to be excluded. In j]^™^"'^''-' vacuum subdiagrams (plain lines) with 4 external lines are 
logarithmically divergent and need to be renormalized. Their sum defines the renormalized four point function A. If 
the latter include the external points (full dots) they contribute to S'"^, defined in (^). All other terms, like the last 
term, contribute to E^"^-*. 



III. THE TADPOLE APPROXIMATION 



The tadpole approximation is given by 

= CO 

.S<^[G] 



3J 



■iT." = 2i- 



5G 



Q 



iA \ - 

A f 



Al_G-(/) 



3J 



(2vr) 



(27r)' 



^16*) 



(17*) 



for the un-renormahzed ^-functional and self-energy. Here a = diag(l, —1) accounts for the inte- 
gration sense on the contour. The self-energy is diagonal in the contour indices. Since the tadpole 
self-energy is real and constant the vacuum part vanishes 



J. (vac) ^ Q 



(18) 



in view of the physical renormalization condition (^) . To find the renormalized equations of motion 
at finite temperature we need the renormalized four-point function A according to (10). Due to 
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(11) we have a constant Bethe-Salpeter kernel 



p(4;vac) ^ ^ ^^^^ 



and thus A is also constant and determined from its renormalization condition (12) 

A^™'^) = ^. (20) 

In the following it is sufficient to just consider the { }-propagator. Writing = + from 

(^ the self-consistent { }-propagator is given by 

^''^P'^ = 2 L ^ ■ - -^^^nripoMp' - m2) 

1 1 1 rj ^ ^ 

with S„{x) = — Im = 7; ^ 

vr X — irj vr + ry^ 

and the Bose-Einstein function 

= SM^i (22) 
Since the tadpole self-energy is constant it is identical with T,^^^ and due to (0) we thus have 
G~W(p) = G"~(p) - G"-(™^)(p) - e[G"-(™")(p)]2 



0^ „ . , ..,2 ,.2. (23) 



(m^ — p"^ — b])^{M^ — — ir]) 

This analytically given result explicitly illustrates that the subtracted part G ^^'^ of the propagator 
no longer obeys the finite temperature KMS and retarded relations (Eqs. (A.14)-(A.18) in I). 
The renormalized effective mass follows from ( [To| ) and (^) with standard integrals from vacuum- 
perturbation theory (see for instance flTl] ): 



= + e = + 



327r2 



x r,,,.!M^^^ ,24) 

ir^ Jo -J I? + hp 



The solution of this gap-equation for m = 0.14 GeV and A = 30 is shown in Fig. ||. The self- 
consistent treatment reduces the mass gap relative to the perturbative result due to the larger 
mass M entering the self-consistent loop. 

The gap equation ( p4[ ) becomes singular for m — > due to the chosen on-shell renormalization 
condition (^. For the sake of completeness we give a brief summary about the treatment of the 
case of vanishing renormalized vacuum mass in the following. 

For m = the renormalization description has to be changed to a so-called "mass independent" 
renormalization scheme. This concept is most easily established in the dimensional regularization 
procedure by the so called minimal subtraction (MS) or modified minimal subtraction (MS) schemes 
where for all renormalization parts, i.e., in our case of (/)'^-theory the proper n-point vertex functions 
with n < 4, only the singular terms in e = (4 — (i)/2 are canceled. 

In the present case of the tadpole approximation where the self-energies are momentum- 
independent this can be done analytically. Indeed for m = we have = and the divergent 
part of the diagram reads 
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Effective Mass for ^=30 (Tadpole) 



> 




0.05 0.1 0.15 0.2 0.25 

T[GeV] 

FIG. 1: The solution of the gap equation ( p4[ ) compared to the naive, i.e., non-resummed one-loop pertur- 
bative solution. Both calculations are renormalized in the vacuum in the same way according to the on-shell 
conditions ^ with m = 0.14 GeV 



in dimensional regularization. Since now the renormalized vacuum propagator is the free massless 
propagator, i.e., 

1 



Gvac(0 



(26) 



+ ir] 

an infrared regulating mass scale fi' has to be introduced in order to calculate both, the mass and 
the four-point vertex 



©vac 

4i(o) 



iA/_f 



2e 



2 



1 



/2 



d'^l 1 



32^^ 



2 J (27r)^ (/2 - ^'2 + i^)2 
with the MS-counter terms 



327r2 



+ 1 - 7 - In 



7 — In 



12 



47r/i2 

12 



+ 0(6) 



+ 0{e) 



6m 



6X 



lim 



A^ 



f2 



327r2 



1 

- -7 



0, 



A 



327r2 



1 



7 



(28*) 



which are IR-finite due to fact that in dimensional regularization the mass counter term in the 
MS scheme is proportional to the infrared-regulator mass /i'^. For m = the renormalized gap 
equation finally takes the form 

AM2 



327r2 

implying that perturbatively 



1 - In 



/r2 



M, 



M 
47r^2 



A 



+ 



A 
4^ 



dLL 



,nr(VL2 + M2) 
VL'^ + M2 



pert 



dLLnriL) 



24' 



(29) 



(30) 



which follows from ( |29D by letting M — > on the r.h.s. 
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IV. THE SUNSET APPROXIMATION 



In this section we calculate the self-consistent self-energy for the next-to-leading order appro- 
ximation of the ^-functional numerically. First we have to solve the renormalized Dyson equation 
for the vacuum and the Bethe-Salpeter ladder equation for A^™'^) , which is needed as input for the 
temperature dependent calculation according to (|lO|). 

The <^-functional is given diagrammatically by 




while the corresponding self-energy and BS-kernel become 



Q 



+ (32*) 

r^') = X + XZ> • (33*) 

For the renormalization parts it is sufficient to restrict all considerations to one real-time con- 
tour branch. For definiteness we choose the time-ordered branch, i.e., in the following subsections 
all contour two-point functions denote { }-quantities rather than matrices. The physical renor- 
malization scheme (^) implies that the tadpole part of the self-energy is already subtracted such 
that in the vacuum case we only need to solve for the self-consistent sunset self-energy. 

The general strategy will be to combine the BPHZ-renormalization scheme with dimensional 
regularization and use the spectral representation for the Green's functions (Lehmann representa- 
tion) 

rd(^ImG'-'H^)^ (34) 

Jq vr — — iry 

All quantities like the renormalized vacuum self-energy J]^™'^^ the BS-kernel F^^^ and the four-point 
function A are then to be expressed through the vacuum spectral function p = —2ImG^^^^^ and 
corresponding time-ordered kernel functions Ki. The latter express the analytic structure entirely 
in terms of free particle properties, however with varying masses. They can then be renormalized 
by dimensional regularization. 

A. The two-propagator loop 

A central quantity is the simple loop-function contained in ( |32*| ) and ( |33*| ) 

(35) 



L(-*5)(g2) = i>Q< =i I _^G(™^)[(/ + <^)2]G(-)(^2)_ 



This four-point function is logarithmitically divergent and to be renormalized with the condition 

L(''^>^)(0) =0. (36) 
With help of the Lehmann representation (B3) it can be expressed as 



^{ren)(^2)^/ ^ ^xf'''^) (g^^ m?, m^) Im G^™^) (m?) Im G^^^^) (m^) (37) 
Jo Jq 
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with the renormahzed kernel defined through 

(38) 

{2'kY — — \rj) [ml — (/ + p)^ — i??] ' 

Here the standard notation for dimensionahy regularized quantities is used, where d = 4 — 2e is 
the space-time dimension, and fi denotes the renormalization scale. 

After a Feynman parameterization the integral (^8|) can be expressed in terms of standard 
formulae of dimensional regularization (see, e.g., |l^]) which leads to the result 




7 s + 



artanh 



s + mf 



m. 



+ {m\ — mi) In 



1712 J 



A(s, mf, mi 
mim2 




artanh 



s — ml + mi 
A(s,mf,mi) 



+ 



with the Kallen function 



A(s,mf,mi) = \J s"^ + mf + m^ — 2sm| H — 2smi — 2m^m2, 



(39) 



(40) 



which is completely symmetric in its 3 arguments and the branch is determined by the s + iry pre- 
scription. Here and in the following we neglect contributions of order 0(e). For the renormalization 
of this kernel we also need its value at s = given by 



Kl'''^\0,mi,m2) 



16(m^ — mi)7r2 



mt 



2 



1 + 7 



+ m^ In 



m^ 



m\ In 



mn 



(41) 



According to (|3q ) we thus find 
<°°)(s,mi,m2) = [<°^^(s) - Kf'^*^) 




mlf 



s{ml + mi) 



m 



In — 



mi 
m2 



(42) 



artanh 



s + mf 



mi 



A(s, m^, mi) 



+ artanh 



s — m\ + mi 
A(s, m\, mi) 



The result for the loop L[s) with s = \s given in Fig. ||. The imaginary part is essentially 
determined by the two-body phase-space which opens at s = 4m^. To discuss the threshold 
singularity of L it is sufficient to study the perturbative result obtained from the pole term of 
Q^^^) which is analytically given as 



L^'^^''\s,m^)=Kr''\s,m\m^) = 



4m ^ 



1 

8^ 



4m2 



■ artanh 



In 



+ V s — 4m2 



2m 



(43) 
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Re L 



Im L 





DC 




9 10 



FIG. 2: Real (left) and imaginary part (right) of the loop function L. For the real part the dashed line 
shows the approximate behavior (^) around the threshold. 

This function is analytic and real for real s < 0. The analytic continuation to other values of s is 
given by taking the principle branch of the square root and the artanh functions for s + irj with a 
small r] > 0. Close to threshold s ~ Am? one obtains the following approximate forms 



ivr Is — Aim? Am? ^ f f s — Am^ \'^\ ^ , r. 

i YV^ r + )) ^<s-^rrP«s. 



(44) 



Around the threshold the square root term causes the singular behavior changing from real to 
imaginary values, while the remaining log-term is entirely real. At large s the real part of L 
behaves logarithmically, whereas the imaginary part goes to a constant. 

B. The vacuum self-energy 



The un-renormalized expression of the two-loop self-energy reads 



E(-=)(p^) = -\j^aj ^G(-^)(/?)G(-=)[(/i + h +pf]G^-\ll). (45*) 



With help of ([37|) the sunset self-energy with all sub-divergences subtracted becomes 



S(™^)(p) = -- 



iA2 r d% 



:L[{h+pf]G^''^^\ll). 



6 J {2ttY^ 

Now we apply the spectral representation for (|35[) with one subtraction determined by 



d(m3)^ 

Am? 



Im L{m'^ + irj) 



1 



1 



irj mg — 11] 



(46) 



(47) 



Together with the Lehman representation of the propagator ( |34D this leads to the renormalized 
vacuum self-energy 
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where due to the renormahzation conditions (^) the kernel ^2'^^'^^ is given by 



K. 



(ren) 



--K 



(rcg) , 



5,^3, 7714) 



K 



(rcg) 



/ 2 2 2\ 



(49) 



The cancellation of the contributions from the subtraction of the sub-divergences is due to a spe- 
cialty of the sunset-diagram: Here all contracted diagrams are of tadpole structure and therefore 
independent of the external momentum p and thus are completely canceled by the overall subtrac- 
tions. The advantage to take them nonetheless into account is that at any stage of the calculation 
we use renormalized functions which can be calculated numerically without using any intermediate 
regularized functions. 

Prom (p9|) we find the analytical expression for the kernel K2: 



K. 



(ren) 



1 



+ artanh 



w?X{s, 7773, 7^-4) X 



7773 -|- m\ 



artanh 



s + mi 



X{s,m3,m4) 



sX{m'^ ,m3, 7774) X 



artanh 



+ artanh 



777 



7n,| + ml 



A(7772, 7773, 7774 




777^ -|- 777^ 



777^ 



X{m'^, 7713, 1^4) 



(50) 



+ 



{s — 777^)[(7n,| — 777|)^ — m^ {m^ + ml)] 



167r'^m^X{m'^,m3, 7714) 



artanh 



771^ + 777I 



777^ 



A(7n^, 7773, 7774) 



+ artanh 



777 



7773 + 7774 



X{m'^ ,m3,m4^) 



777 



2\2 



+ 



(7773 — 7774)(s — 777 ) 



■In 



The equations (p^) and (I 



Idn'^m'^ IGvr^ 
supplemented by the Dyson equation 



7773 
7774 



1 



p. 



m^ 



S(™=)(p2) +ir,' 



(51) 



form the closed set of renormalized equations of motion, which were solved with the analytically 



given Kernels K^'^^^ and K)^^''' from (|4^) and (|50|). For the integrals (|37| ) and (|48| ) a simple adaptive 
Simpson integrator was used. We have chosen 777 = m-,^ = 140 MeV for the mass and A = 30. As 
Fig. ^ shows for this coupling there is no visible difference between the perturbative and the self- 
consistent result since the main contribution comes from the pole term of the propagator, while 
the continuous part, which starts at a threshold of s = 9777^, is suppressed. Thus, due to our 
renormahzation scheme, where 777 is the physical mass parameter, the pole-term result essentially 
coincides with the perturbative one. 



,(rcn) 



C. The Bethe-Salpeter function 

The renormalized Bethe-Salpeter kernel F^^''''^") defined in ( [HI) reads 

F(^''-^'^)(/,p) := f(^''-^'^)[(Z = ^ + X^L'^'''''\{l-pf]. 



(52) 
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FIG. 3: Real (left) and imaginary part (right) of the sunset self- energy. The perturbative and the self- 
consistent result lie on top of each other due to the large threshold at s = 9to^. 



Since as shown above the pole term essentially determines the loop-function L, we simplified the 
task to solve (|l^) by using the free propagator 



1 



p2 _ jj^I _|_ jj^ 



(53) 



instead of the self-consistent vacuum 

Subtracted dispersion relations for r^^''''^'^) and A^'''^") and the renormalization conditions 
p(4,ren)^Q^ _ A^'''^'^) (0, 0) = A/2 provide the renormalized integral equation 



d(m2) 



4m2 ^ 



Imr(^'''^")(mi) 



2lp - 



{I — — i??](?7T-| — /2 — \rj) 



(54) 



d(m?)/2imA('-™)(0,mf) A 



m 



1 — P —irj 



+ 



Again the ^-integration can be performed with the help of standard perturbation theory integrals 
leading to the kernels 

^2(2/p-p2^ 



r^(ren) / 2 n ■ 

(P ,'"T'l,'"T'2) = 1 

r^(ren) / 2 \ ■ 



dH 



(27r)^ (m^ — P — irjy (mf — P — ir?)[rfi2 — {I — pY — if]] ' 



2lp — p 



(55) 



(27r)4 {m? -P- ir/)2 [ml -{I- pf - ir;] 
which can be related to the previous kernel Ki through the identity 



1 



1 



2m mP —P — if] 

We start with the calculation of defined in ( |55| ) and write it in the form 



2m 



-dr, 



dH 



1 



(27r)^ m^ — P — 17] 



1 



ml — {I — p)^ — nT-l ~ 



17] 



(56) 



(57) 
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FIG. 4: Real (left) and imaginary part (right) of the A-function. The plot shows that after three iterations 
the solution is already stable. 



which can be expressed in terms of 

1 



(ron) 



K4{s, m, 7712) 



2m 



dmKi{s,m,m2) 



+ 



{iv? — m\ — s) 



IGvr^s I — 771^ \{s,m,m2 



. , s + 777? — rrin \ 1 f ^ ~ ""^^ + ^2 

artann ( — ^ I + artann ' 



A(s, 7771, 7772) 



A(S, 7771, 7772) 



777"^ — 1712) + 2m2S 



In 



777 



(777^ — 7773)^ \'m2. 

The integrand of cf. (^) can be rewritten as follows 

i /2 p2 _ 



^ — P — IT] mf (mf — P — ir]) [m'2 — P — irj) [m'2 — {I — p)'^ — i?/] 



777 



1 



-—dm 

2777 777^ — 777^ ^777^ — P — IT] 171^ — P — 17] 



m\{rn? — P — bj) 
The integration of (p^) then leads to 



1712 — {I — P)^ — i?? 7X12 ~ l"^ ~ i^. 



7^^3(8,777, 7771,7772) 



Kl(s, 7771, 7772) — i^l(s, 777, 7772) 777^ — 2777i 



+ 



(777^ — 777^)2 (777^ — m'\)m\ 

With the so defined kernels we can express (p3) as follows: 



A(-'^)(0,p2) = r(p2)+ / ^H!!!2 7^3(^2, 7771, 7772) ImA('-<=°)(0, 777?) Imr(4''-°°)(7n2) + 

74^2 7r777f TT 



A 



d777n 



2 J im? 



K4(p^r772)Im^(^■'"™)(777^) 



(58) 



(59) 



i4r4(s, 777, 7772). (60) 



(61) 



Fig. 1^ shows the solution of this renormalized equation of motion in comparison with the one-loop 
approximation 

p(4,rcn)^ As the fi gure shows, the solution is stable after three iterations. 
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D. The temperature dependent parts 

To calculate the temperature dependent part of the self-energy we use the analytical properties 
of the two-point functions and the fact that the retarded propagator fulfills the simple Dyson 
equation 

GdP) = -2 2^wT' (62) 

where the retarded self-energy can be obtained from the { }-matrix element via the equilibrium 

expression 

Sfi(p) = ReS— (p) +itanh (^) ImS— (p). (63) 

As explained in I both, G^™'^^ and are diagonal matrices while the corresponding off-diagonal 
parts are contained in G'^^\ because they are of a lower degree of divergence due to explicit Q -\-n- 
f actors. 

Due to its topology the sunset-diagram has no further vertices besides its two external points. 
Therefore the 2PI vacuum four-point function F^^'™'^) is effectively only a two-point function and 

thus we need only the { }-components of the various Green's functions in order to calculate 

E . This implies 

Q (matter) q q (vac) ^ (r) q (matter) (vac)j2j-i (0) (64) 

With help of the renormalized vacuum parts we obtain for the self-energy 

s-(p) =s-(-^)(p) + 



+ 1 



/1 4 7 
_£_ |r(4,vac)j(^ ^ ^)2] _ p(4,vac)^^2)| ^—(matter) ^ 



+ ^/(^^'^'^^^°'^')^"'^^^)- (65) 
"y/ (04 / (^G~(")(/i)G~(")(/i + Z2)G-(™^)(/2+p)- 
"y/ Wf I (^G--(--"<''-)(/i)G--(--"°'-)(/i + Z2)G--(--"^'-)(/2+P^ 

Note that all vacuum quantities entering here are to be taken in their renormalized version and 
that the remaining integrals are all finite due to power counting. 



E. Results 

The above finite integrals are to be evaluated numerically. While due to Lorentz invariance 
the vacuum loops involve just two-dimensional integrations which can numerically be integrated 
by standard methods, for the in-matter loops only rotational symmetry in three momentum space 
can be exploited. For each loop diagram this leads to three dimensional integrals for each external 
momentum PO) We solve these integrals on an equal distant N x N lattice in these coordinates. 
Naively the computing effort would then scale like N^. However, we succeeded to develop an 
improved algorithm for the loop integrations, where the computing time scales with N'^ (essentially 
a gain of more than two orders of magnitude). The lattice implies an infrared cut-off which 
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Pert. ReE for T=250MeV, >.=30 Pert. ImE for T=250MeV, 1=30 




ReE for T=250MeV, ^1=30 ImE for T=250MeV, 1=30 




FIG. 5: Real (left) and imaginary part (right) of the perturbative (top) and the self-consistent self-energy 
for A = 30, m = 140MeV and T = 250MeV. Note that the self-energies are multiplied with an factor 100 in 
these plots! 

requires particular care for the treatment of the sharp structures of the vacuum propagator around 
its on-shell pole. The method at hand was to tabulate the lattice cell integrated values of the 
propagator and its moment in pQ given by the analytic result from a linear interpolation in Pq 
of the inverse propagator at fixed \p\. The remaining self-energy factors of the integrands are 
approximated to linear order in Pq. The interpolation procedure is adapted to the fact that due 
to the quadratic pole of the {G (vac)-)2 ^gj,j^ jg sensitive to the po derivative (residue) of the 
remaining integrand. Both, the infrared and the UV cut-offs of the lattice have been controlled by 
varying the lattice spacing and the overall size of the lattice. The final results were achieved with 
Apo = A|p1 = 6.67 MeV with = 300, i.e., a UV-cutoff at 2 GeV. 

Perturbative results are obtained through Eq. ( [65| ) approximating q("^^^^'^'^) by the KMS- 
temperature part of the free propagator G^™'^'^) (last term in Eq. (|l3|); note that here G^''^ = 
^(matter) gjj^gg ^(vac) _ q-j_ The self-cousistent solutions are then obtained iteratively through 
the set of Eqs. (|6^) to (|65|). The results for both, the perturbative and the self-consistent case, 
are shown in Fig. ^ in a 3-dimensional plot over the [po, |p|)-plane illustrating that the entire 
calculations are performed with the full dependence on energy and momentum. Details can be 
extracted from the cuts shown for a set of selected momenta in Fig. ^. The main qualitative 
results are similar for both the perturbative and the self-consistent calculation: In the vacuum and 
self-consistent pure tadpole case the self-energy shows a threshold cut resulting from the decay into 
three particles, i.e., Pq— j?' > 9M^. Adding the sunset self-energy leads to a spectral width which 
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FIG. 6: Real (left) and imaginary part (right) of the perturbative (top) and the self-consistent self-energy 
for A = 30, m = 140MeV and T — 250 MeV as a function of po for various 3-monienta. Note that the 
self-energies are multiplied by a factor of 100 in these plots! 



dissolved this threshold such that the self-energy shows spectral strength (imaginary parts) at all 
energies. While the growing high-energy tail is related to the decay of virtual bosons into three 
particles, at finite temperature, as a new component, a low-energy plateau in ImS^ emerges from 
in-medium scattering processes. 

Various balancing effects are encountered for the self-consistent case: For sufficiently large cou- 
plings and/or temperatures the self-consistent treatment shows quantitative effects on the width. 
The finite spectral width itself leads to a further broadening of the width and a smoothing of the 
structures as a function of energy. This is however counter balanced by the behavior of the real 
part of the self-energy, which, as discussed below, essentially shifts the in-medium mass upwards. 
This reduces the available phase space for real processes. With increasing coupling strength A a 
nearly linear behavior of Im with po results implying a nearly constant damping width given 
by — Im S^/po- 

The overall normalization of the real part of is determined by the renormalization procedure. 
In this case there are three counter balancing effects. First the tadpole loop shifts the mass to 
higher values. As the tadpole is less effective for higher masses this effect weakens itself in the 
self-consistent tadpole treatment, c.f. Figs. |l| or 0. However, since the sunset part adds spectral 
width, it indirectly contributes to the tadpole loop. Since spectral strength at the lower mass side 
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FIG. 7: The in-medium effective masses M (left) and spectral widths T (right) of the particles for the various 
approximations described in the text as a function of the system's temperature T. 



carries higher statistical weights, the tadpole loop in turn leads to a further increase of the mass 
shift, c.f. the perturbative calculations of sunset & tadpole in Fig. |^ 

The direct contributions of the sunset terms to the real part of the self-energy become relevant 
at higher couplings and temperatures. Then the self-consistency leads to significant effects which 
contributes to a net downshift of the real part of the self energy or in-medium mass M. The latter 
effect finally overrules the tadpole shift and indeed leads to an overall negative mass shift compared 
to the (tadpole dominated) perturbative result. These effects are illustrated in Fig. ^ where the 
in-medium effective mass M and width F of the corresponding "quasi-particles" are plotted against 
the temperature. Thereby M and F are defined as the quasi-particle energy M = pQ at vanishing 
real part of the dispersion relation [pg — ttt? — ReT,{po,p)]^^^^,j = for p = 0, and through 

F = -ImS(p)/po|p=(M,o)' respectively. 



V. CONCLUSIONS AND OUTLOOK 

We have shown that it is possible to use the renormalization scheme, proposed in [Q], for 
numerical investigations of the self-consistent approximations for the self-energy derived from the 
truncated effective action formalism on the 2PI level. Thereby it is very important to isolate the 
divergent vacuum parts consistently, in particular the implicit or hidden ones, from all convergent 
and in particular temperature or more generally matter-dependent parts. This could be provided 
by the ansatz given in I for both, the propagator and the self-energy. The renormalized vacuum 
pieces are obtained using the Lehmann representation for all two-point quantities. The resulting 
integration kernels can then be renormalized by standard technics. The procedure solely rests 
on Weinberg's power counting theorem, i.e., on an analysis of the asymptotic behavior of the 
propagators. 

In this way both, the renormalized vacuum self-energy and four-point functions can directly be 
obtained from finite equations. For the finite temperature parts it is important that they have 
to be completely excluded from the counter term structures. This is achieved by the technique 
developed in I. Exploratory calculations are shown for the symmetric Wigner-Weyl phase of the 
X(p^ model for the self-consistent treatment of both the tadpole and the sunset diagram at finite 
temperature. 

The results promise that the method, which is conserving M, 03] and thermodynamically con- 



sistent, can also be applied for the genuine non-equilibrium case, i.e., in quantum transport |18] 
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or for the solution of the renormalized Kadanoff-Baym equations. Numerical studies for non- 
equilibrium cases (Gaussian initial conditions, spatially homogeneous systems) were already per- 
formed in [|l^, |2^, These investigations were undertaken in 1-1-1 dimensions. Our renor- 
malization scheme should also be applicable for the non-equilibrium case in 1-1-3 dimensions and 
implement able in numerical codes. 

The investigation of the symmetry properties of <I>-derivable approximations is the subject 
of a forthcoming publication 0. It is known that in general the symmetries of the classical 
action which lead to Ward-Takahashi identities for the proper vertex-functions are violated for 
the self-consistent Dyson resummation for the functions beyond the one-point level, i.e., on the 
correlator level. The reason is that, although the functional T can be expanded with respect to 
expansion parameters like the coupling or h (loop expansion) or large- expansions for O(A^) type 
models, the solution of the self-consistent equations of motion contains partial contributions to 
any order of the expansion parameter. This resummation is of course incomplete and violates 
even crossing symmetry for the vertices involved in the renormalization procedure. This causes 
problems concerning the Nambu-Goldstone modes Q in the broken symmetry case or concerning 
local symmetries (gauge symmetries) 1^] on a level where the gauge fields are treated beyond the 
classical field level, i.e., on the propagator level. 

It can be shown though, that on top of any solution of a ^-derivable approximation which is 
constructed from a symmetric Lagrangian there exists a non-perturbative effective action regf[y5] 
which generates proper vertex functions in the same sense as the IPI effective action. These ex- 
ternal vertex functions fulfill the Ward-Takahashi identities of the underlying symmetry. However, 
in general they coincide with the self-consistent ones only up to one-point order. This fact espe- 
cially ensures that the expectation values of Noether currents are conserved for the ^-derivable 
approximations. Thus usually the so generated external self-energy and higher vertex functions 
are different from the <I>-derivable expressions. Details on these considerations will be given in a 
forthcoming paper 
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